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INTRODUCTION 


The  Air  Force  has  a  vital  interest  protecting  sensitive  electronic 
equipment  inside  airplanes  from  electromagnetic  pulse  (EMP)  penetration 
through  composite  skin  panels.  Composite  skin  panels  are  being  used  on 
airplanes  today  instead  of  aluminum  because  they  are  stronger  and  lighter 
than  aluminum.  Composite  panels  offer  less  shielding  than  conventional 
metal  panels,  however.  Besides  direct  penetration  (diffusion)  through 
the  composite  panel  walls,  inadvertent  penetration  through  seams ,  joints, 
and  windows  is  possible. 

The  evaluation  of  the  electrical  (shielding)  properties  of  composite 
materials  has  been  studied  by  Casey  (1976,  1977).  For  a  graphite  composite, 
he  concludes  that  the  electrical  conductivity  in  the  direction  normal  to 
the  sheet  surface  is  approximately  zero.  The  tangential  conductivity  is 
anisotropic,  but  he  shows  that  it  can  be  satisfactorily  modelled  by  an 
isotropic  conductivity  of  approximately  1 . 5xl04mhos/meter  (for  graphite). 

He  further  shows  that  the  graphite  composite  acts,  in  effect,  like  a  low 
pass  filter.  This  means  that  only  lower  frequencies  usually  need  to  be 
considered. 

Because  the  frequencies  of  interest  are  low,  a  sheet  of  graphite 

composite  material  of  thickness  6  and  effective  tangential  conductivity 

a  can  be  modelled  by  an  infinitely  thin  sheet  of  sheet  impedance  Z  =  l/o  5. 
t  sc 

Ohm's  law  requries  that  inside  the  sheet, 


where  E  is  the  tangential  component  of  electric  field  in  Volts/meter 
^  tan 

and  J  is  the  sheet  current  in  Amps  per  meter. 

This  report  describes  a  method  of  solving  for  the  deterioration  of 
shielding  due  to  slits  (seams)  in  composite  materials.  This  is  done  by 
solving  for  the  quasi-static  (low  frequency)  magnetic  field  diffracted 
by  a  slit  in  an  impedance  sheet.  Only  the  E-polarization  is  studied.  Three 
different  electromagnetic  concepts  are  involved  in  treating  this  problem. 
These  are 


(a)  Low  frequency  or  quasi-static  techniques 

(b)  Impedance  boundary  conditions 
and  (c)  Scattering  by  slits. 

Previous  efforts  which  combine  these  three  topics  have  been  few.  Hurd 
(1979)  has  recently  treated  a  similar  problem  for  the  H-polarization.  He 
uses  an  impedance  plane  instead  of  an  impedance  sheet  boundary  condition. 
Impedance  planes  have  a  surface  impedance  boundary  condition  whereas  impedance 
sheets  have  a  jump  discontinuity  boundary  condition.  Kaden  (1959,  p.  212) 
uses  a  conformal  transformation  to  find  the  penetration  through  a  gap  in  a 
plane  shield  with  finite  conductivity  and  thickness.  This  work  has  been 
summarized  by  Butler,  et  al  (1976). 

The  literature  that  describes  one  of  the  three  subjects  individually 
is  reviewed  briefly  below. 

(a)  References  on  low-frequency  techniques. 

Low  frequency  scattering  techniques  are  often  used  because  analytical 
solutions  can  sometimes  be  obtained.  Two  review  articles  on  low  frequency 
techniques  have  been  written  by  Kleinman  (1967,  1978).  Quasi-static 
techniques  reduce  the  dynamic  electromagnetics  problem  to  a  static  problem 
by  making  suitable  low  frequency  approximations.  Latham  and  Lee  (1968) 
develop  quasi-static  boundary  conditions  for  inductive  shields  by  neglecting 
displacement  currents.  Standard  techniques  for  solving  statics  problems 
can  then  be  applied.  Such  techniques  are  detailed  by  Sneddon  (1966),  among 
others. 

(b)  References  on  impedance  boundary  conditions. 

Use  of  impedance  sheet  boundary  conditions  has  been  of  recent  interest. 
Harrington  and  Mautz  (1975)  use  them  to  treat  thin  dielectric  shells.  Senior 
(1978)  discusses  them  in  connection  with  impedance  half  planes  and  Senior 
(1979)  discusses  them  in  conection  with  resistive  strips.  Casey  (1977)  gives 
conditions  under  which  composites  can  be  modelled  by  impedance  sheet  boundary 
conditions.  Babinet's  principle  for  impedance  boundary  conditions  has  been 
given  by  Baum  and  Singaraju  (1974),  Lang  (1973),  and  Senior  (1977). 

(c)  References  on  scattering  by  slits. 
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Scattering  by  slits  in  perfectly  conducting  planes  has  been  studied 
exhaustively.  Scattering  by  slits  in  finitely  conducting  planes  or  in  impedance 
sheets  has  not  received  much  attention.  Slits  in  perfectly  conducting  planes 
have  been  studied  by  Clemmow  (1966),  Hongo  (1972),  Houlberg  (1967),  Nomura  and 
Katsura  (1957),  Millar  (1960),  and  Otsuki  (1976),  among  others.  Lam  (1976) 
has  studied  the  shielding  effectiveness  of  seams  or  joints  in  perfectly  con¬ 
ducting  aircraft  skins.  Diffraction  by  slits  or  apertures  in  impedance  planes 
has  been  treated  by  Neugebauer  (1956),  Zakharyev,  Lemanski  and  Shcheglov  (1970), 
and  Hongo  (1972).  Neugebauer  develops  an  approach  to  treating  apertures  in 
absorbing  screens  by  using  symmetry  conditions.  Zakharyev,  et  al,  handle  the 
case  of  a  dipole  antenna  located  in  a  slit  in  a  finitely  conducting  plane. 
Finally,  Hongo  uses  the  Weber-Schafheitlin  integral  to  formulate  the  problem 
of  diffraction  by  a  slit  in  a  screen  with  a  surface  impedance.  Apparently, 
no  one  has  previously  studied  the  probelm  of  interest  here  —  diffraction  by 
a  slit  in  an  impedance  sheet. 

IT.  OBJECTIVES 

The  objective  of  this  research  is  to  find  the  H-field  diffracted  by  a 
narrow  slit  in  an  impedance  (composite)  sheet.  Equations  which  are  functions 
of  slit  width  and  panel  thickness  and  effective  conductivity  are  to  be 
developed  and  solved.  The  solutions  can  then  be  used  to  study  the  deterioration 
of  shielding  due  to  slits  (seams)  in  composite  panels.  Three  possibilities 
exist  for  each  slit  width  and  panel  thickness  and  conductivity.  One  might 
find  that  penetration  through  the  slit  is  the  major  component  of  coupling. 

On  the  other  hand,  one  might  find  that  penetration  through  the  slit  can  be 
neglected  compared  to  the  direct  penetration  (diffusion)  through  the  walls. 
Finally,  it  might  turn  out  that  both  coupling  through  the  slit  and  the  walls 
has  to  be  considered. 

III.  REFLECTION  AND  TRANSMISSION  COEFFICIENTS  OF  AND  THE  CURRENT  IN  AN 

IMPEDANCE  SHEET 

The  study  of  an  impedance  sheet  without  a  slit  will  be  a  helpful 
start  and  of  use  later.  A  similar  discussion  can  be  found  in  Seshadri 
(1971).  Figure  1  shows  an  E-polarized  plane  wave  incident  upon  an 
impedance  sheet  at  an  angle  9.  The  sheet  is  in  the  x  =  0  plane. 


Figure  1.  Plane  wave  incident  on  impedance  sheet. 
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where  o  /£ o  .  There  are  two  boundary  conditions  that  must 

be  satisfied  across  the  impedance  sheet  [Casey  (1977,  p.  23)].  These  are 
that  E  tangential  (E  )  must  be  continuous  and  that 

Lan 

Eh>-o)  —  2S  J  =  )f) 

Applying  these  conditions,  one  obtains 
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This  current  satisfies  the  integral  equation 
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as  may  be  easily  shown  by  a  Fourier  transformation. 
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IV.  SLIT  INTEGRAL  EQUATIONS:  DYNAMIC  CASE 


A  standard  approach  for  solving  electromagnetics  problems  is  the 

integral  equation  approach.  Consider  a  two-dimensional  slit  of  width 

2b  between  two  impedance  sheets  of  impedances  Z  and  Z  ,  as  shown  in 

S1  s2 

X 


Figure  2.  Dynamic  plane  wave  incident  on  slit  in  impedance  plane. 
Figure  2.  The  vector  potential  for  this  problem  is 


/-CO 

+  if- 

K 

The  electric  field  scattered  by  the  impedance  sheets  is  given  by 

Es  -  (-1. 

f 

Since  the  total  electric  field  is  equal  to  the  incident  field  plus  the 
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scattered  field  and  is  also  equal  to  Z  I  over  <  z  <  -b  and  to  7.  1_ 

SI  1  s2  2 

over  b  <  z  <  °°,  both  for  x=0,  one  obtains  the  coupled  integral  equations 
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These  equations  can  be  solved  by  using  numerical  techniques. 

Since  these  equations  are  very  complicated,  it  is  desirable  to 
simplify  them.  One  possibility  would  be  to  try  using  n  Babinet’s  prin¬ 
ciple  of  the  types  described  by  Baum  and  Singaraju  (1974),  Lang  (1973), 
or  Senior  (1977).  Unfortunately,  this  does  not  lead  to  a  significant 
simplification  of  the  equations  and  so  is  not  considered  further.  Another 
possibility  is  to  look  at  the  equations  for  low  frequencies.  This  also 
does  not  lead  to  a  major  simplification  in  the  equations  as  they  stand 
because  the  currents  are  supported  over  two  semi -inf  ini tc  regions  and 
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and  are  therefore  not  directly  tenable  to  low  frequency  approximation:; . 

A  solution  can  be  found  using  this  method,  however,  by  solving  for  the 
difference  between  the  actual  currents  and  the  currents  that  would  be 
present  without  the  slit.  This  perturbation  current  goes  to  zero  as  one 
gets  away  from  the  slit.  This  approach  is  the  one  taken  here. 

V.  QUASI-STATIC  DUAL  INTEGRAL  EQUATION  FOlCMULATION  OF  THE  SLIT  PROBLEM 
The  problem  that  will  be  studied  here  is  less  general  than  the  one 
detailed  in  the  last  section.  First  of  all,  it  will  be  treated  on  a 
quasi-static  basis.  Secondly,  only  the  case  of  normal  incidence  will  be 
studied.  Finally,  both  parts  of  the  sheet  have  the  same  impedance.  This 
is  shown  in  Figure  3. 


Figure  3.  Quasi-static  field  incident  on  slit  in  impedance  plane. 

-V 

Let  the  incident  field  H  =  H^z  be  the  field  of  a  static  current 
sheet  of  strength  2Hq9  located  at  infinity.  It  is  of  interest  here  to 
find  the  field  above  and  below  the  impedance  sheet  and  also  in  the  gap. 
The  basis  for  the  quasi-static  formulation  used  here  is  provided  by 
Latham  and  Lee  (1968).  Maxwell's  equations  can  be  written  ns 


Cs.o 
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(5.4) 


Thus,  these  equations  can  be  written  as 


V7*  V\  =  J 


V  x'E 


^  CO  fA  |_| 


(5.0) 


provided  that  Zg  is  small  enough  and  the  frequency  is  low  enough  so  that 
the  displacement  current  can  be  neglected.  Off  of  the  sheets,  one  has 


=  o  (5  •“if 

v  -  H  ~  O  (5  ?) 


The  first  condition  implies  that 

H  =  -vtf  .  fs.v) 

where  <J>  is  the  magnetic  scalar  potential  [Stratton  (1941),  pp.  225-2671 
while  the  second  gives  that 

v2  $  =  o  (5.10) 
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away  from  the  sheets.  The  incident  potential  is  related  to  the  incident 
field  by 

H'  =  -V0‘  (5.(1) 

so  that 

=  -  Ho  'Z  ■  (5. 12) 

The  arbitrary  constant  of  integration  is  chosen  to  be  zero  so  that 
<J>1(z=Q)  =  0.  This  makes  the  incident  potential  an  odd  function — a  property 
which  will  be  very  useful  later. 

It  is  now  necessary  to  determine  the  boundary  conditions  which  the 
potential  must  satisfy.  Latham  and  Lee  (1968)  show  that  the  surface  diver- 
gence  of  J  is  zero,  that  is, 

Vs  •  J  =  O  (5.1 1> 

in  the  quasi-static  case  for  an  impedance  sheet.  Thus,  the  current  can  be 
written  as 

J  =  ■y.  «  (5.14) 

where  iji  is  a  scalar  function.  This  follows  from  equation  38,  page  502  of 
Van  Bladel  (1964).  Figure  4  shows  the  cross-secional  view  of  the  imped¬ 
ance  sheet  and  a  surface  S  of  integration.  Integrating  equation  (5.5) 


44 


Figure  4.  Surface  and  path  of  integration, 
over  the  surface  S  gives  » 
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=  J>bi-dL  =  p-  H-,(o+)  +  f-lj(cr) 


A2 


=  1  J'dS  =  (- (5.15) 

S 


After  using  (5.9),  one  obtains 


0(0+) 


0(0")  =  -'Y.  (5.10 

^  H  •  = 


This  relates  the  potential  to  i{».  Since  /A  OP  H  *  d  —  O  over  a  small 
"pillbox"  enclosing  the  surface,  normal  H  is  continuous  across  the  surface 


and 


<3$ 

3$ 

^x 

1+ 

31 

3x 


(5  n) 


on  the  boundary.  One  last  condition  must  be  found.  Latham  and  Lee  (1968, 
p.  1749)  show  that 

a.<jo/a  l-U  —  %  •  Vs  *  E^  =  *£  *  X7S  x(?sT^ 


so  that 


-AcojAr-'xnfi  =  Zs  ^  *  (ys  *  ^  (e.  I  ?) 


Therefore,  one  has 


26 

ACOyfA  ^-x 


=■  Zs  (^z  t) 


(S.vcA 
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after  using  identity  37,  page  502  of  Van  Bladel  (1964).  Eliminating  l|> 
between  this  equation  and  (5.16),  one  arrives  nt  the  second  boundary  con¬ 
dition  that  must  be  satisfied,  namely, 


(COM  HI  (\7^(  0(of,2)  -  7$))  (5.21) 

o 

One  might  expect  that  a  canonical  problem  to  test  these  boundary 

conditions  would  be  that  of  scattering  by  an  impedance  sheet  without  a 

slit.  This  does  not  lead  to  any  useful  result,  however.  Instead,  one 

can  use  the  results  of  section  III  for  normal  incidence.  Thus,  for  k**0 

and  G  =  90°,  one  has  that  H^HgZ,  HC  =  THQz ,  and  Hr  =  -R'HQz  where  T  - 

2Z  /(Z„+2Z  )  and  R'  =  -Zn/(Z„+2Z  ).  It  can  thus  be  assumed  that  for  = 
s  0  s  0  0s 

-v 

=  -THoZ  .m 

where  R  =  -R’.  By  using  these  reflected  and  transmitted  potentials,  one 

g 

can  formulate  the  problem  in  terms  of  a  potential  <f>  which  will  be  signif¬ 
icant  only  in  the  neighborhood  of  the  slit.  For  convenience,  define 


This  potential  must  satisfy  Laplace's  equation  (5.10).  Using  separation  of 
variables  techniques  and  integrating  over  all  possible  separation  constants. 


(s.V'M 
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The  total  potential  can  be  written  as 


0 


-H 


where 


and 


0"  =  0C  1-  (f> 

Applying  the  condition  (5.17),  one  obtains 


F  M  =  *  F2  (x) 


(5.25) 

T<o 

(s.  Zfco) 

(5.2U) 

:ains 

)  =  -  f(\) . 

(5-.  27) 

Substituting  this,  (5.24),  (5.22)  in  (5.26),  one  arrives  at 


0+  =  ^  e  V*  4  -  FoZ  C-5.7^) 

0-  =.  -K01  FO-HivA*  e'W4  -  T  Mo?  C5.-ZS0 

■'o  ^ 


It  is  now  necessary  to  find  F(X) .  To  do  this,  the  previously  derived 
boundary  condition  (5.21)  can  be  applied  on  the  impedance  sheet.  In  the 
slit,  the  magnetic  scalar  potential  must  be  continuous.  The  application 
of  these  two  conditions  gives  dual  integral  equations  which  can  be  solved 
for  F(A).  Application  of  the  first  condition  gives 

cO 


A 

Aiz 


-MO/* 


9) 
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,  we  have 


2Z, 


Since  1+R-T  = 


Z-+2Z 
0  s 


-  0~(o,l)  = 


Jk 

T 


S'UU 

2)-l0\  FO)«**>Xz 

"2^o  ? (20+227) 


(s-3o) 


and  (5.29)  becomes 


roo 


( S.Vc *) 

1?  I  >  b 


Application  of  the  second  condition,  continuity  in  the  slit,  yields 


(s.30) 


These  two  equations  are  the  dual  integral  equations  that  must  be  solved  for  the 
unknown  F(A).  Once  F(  A  )  is  known,  the  potential  can  be  determined  from 
(5.28). 

The  following  section  discusses  the  evaluation  of  some  integrals  that 
we  will  need  to  solve  the  problem  and  the  section  after  that  solves  the 
equations  in  the  case  of  a  perfectly  conducting  screen.  Section  VIII  gives 
several  possible  formulations  of  the  dual  integral  equations. 


VI.  generalized  HANKEL  TRANSFORMS 

In  the  process  of  formulating  and  solving  this  problem,  one  finds  it 
necessary  to  evaluate  a  great  many  integrals.  Many  of  these  integrals  can  be 
expressed  as  conventional  Hankel  transforms.  If  one  denotes  the  class  of  all 
absolutely  integrable  functions  f(x)  on  the  interval  (0,  OO  )  by  L^(0,  00  ) , 
each  f(x)  £  1^(0,  00  )  must  satisfy 
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jjlMlc/vr  <<£>  .  (£-  ^ 

If  f (x)  €  Lj(0,  00  ).  then  the  conventional  Hankel  transform  is  given  by 

F(X)  =  (fc.Z) 

where  (x)  is  the  Bessel  function  of  the  first  kind.  If  an  integral  can 
be  cast  in  this  form,  it  is  easily  looked  up  since  several  integral  tables 
of  conventional  Hankel  transforms  are  available. 

In  the  following,  it  will  occasionally  be  necessary  to  evaluate  the 
Hankel  transform  integral  when  f(x)  is  not  absolutely  integrable.  If  one 
attempts  to  look  such  an  integral  up  in  a  table,  they  will  often  find  that  the 
limits  of  applicability  given  next  to  the  formula  are  violated.  In  other 
cases,  one  will  find  that  the  value  of  the  integral  is  listed  as  infinity. 

This  can  be  particularly  disturbing  If  the  integral  was  supposed  to  be  a 
kernel  or  a  right-hand  side  in  an  integral  equation. 

The  solution  of  the  slit  problem  at  hand  is  based  in  part  on  the  classical 
theory  of  dual  integral  equations  presented  by  Sneddon  (1966) .  This  book 
is  based  on  the  conventional  Hankel  transform  as  is  a  later  publication  of 
his  given  in  Ross  (1975).  The  use  of  the  generalized  Hankel  transform  to 
formulate  and  solve  dual  integral  equations  has  apparently  not  received  a 
great  deal  of  attention  to  date  except  by  Zemanian  (1968)  and  Walton  (1973). 

The  following  discussion  of  the  generalized  Hankel  transform  is  essentially 
a  brief  summary  of  that  given  by  Zemanian  (1968) . 

There  are  two  equivalent  ways  in  which  a  generalized  function  may  be 
defined.  First,  a  generalized  function  may  be  defined  to  be  a  continuous  linear 
functional  on  some  space  of  testing  functions.  Second,  a  generalized  function 
by  be  defined  as  a  limit  of  equivalent  fundamental  sequences  of  continuous 
functions.  To  each  generalized  function  in  the  functional  approach,  there  is 
one  corresponding  generalized  function  in  the  sequential  approach.  In  order  to 
follow  Zemanian  (1968),  the  functional  approach  is  used  here. 

In  order  to  define  a  functional,  one  must  first  define  a  space  of 
testing  functions.  For  the  generalized  Hankel  transform,  Zemanian  chooses 
a  class  of  testing  functions  ^  of  rapid  descent.  A  function  too  is  in 
if  and  only  if 
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a.  it  is  defined  on  0  <  x<rco 

b.  it  is  complex  valued  and  smooth 

c.  all  derivatives  are  of  rapid  descent. 

For  the  exact  mathematical  meaning  of  these  terms,  the  reader  is  referred  to 
Zemanian. 

DEFINITION :  A  generalized  function  is  a  continuous  linear  functional  f  on  a 
testing  function  space  (in  this  case  %>■ 

A  functional  f  on  ^j^is  a  rule  that  assigns  a  complex  number  to  every 

member  .  A  continuous  linear  functional  f  on  %  possesses 

(1)  linearity;  that  is,  for  any  two  testing  functions  ^  and  *n 

and  any  complex  number  cK 

<i,  0,+^?  =  <r-P,  a, y  +<r-f, 

<*  <fi~>  ~  y 

(2)  continuity;  that  is,  for  any  sequence  {  (Jr)  ^  in  n- 

that  converges  in  to  ft 

f  tv]^  ^  • 


f  *  in  U*  • 


The  collection  of  all  continuous  linear  functionals  on^^,  is  called  the  dual 
space  of  %  and  is  denoted  by  •  The  generalized  functions  f  are 
sometimes  called  generalized  functions  of  slow  growth. 

A  locally  integrable  function  f(x)  is  said  to  be  a  function  of  slow 
growth  if  a  .  .  ,-N 


JL*  1 1 1  N  -f  (4)  -  0 


(&y 


for  some  integer  N.  Let  f(x)  be  a  locally  integrable  function  on  0  <£x<QD 
such  that  f(x)  is  of  slow  growth  as  x-yoo  and  x^+ls  f(x)  is  absolutely 
integrable  on  0  <  x  < 1.  Then,  f(x)  generates  a  regular  generalized  function 
f  in  *W^by 


J  |(v;  too  b 


fc-4) 


The  generalized  Hankel  transform  is  based  on  Parseval's  formula  for 
conventional  Hankel  transforms.  If  f(x)  and  G(y)  are  in  L^(0,  00  >,  if  //  z-£  , 
if  F(y)  -  and  if  g(x)  =  K'  ,  then  the  Parseval  formula  states 


is  the  symbol  for  the  Hankel  transform.  If  we  let  G  be  an  arbitrary 


$  €  such  that  and  ^et  f  be  an  arbitrary  f  £  ’  we  ^effne 

the  generalized  Hankel  transform  F  =  f  by 

<F/  §>  >  -  <’-f/ 

fa) 

By  using  this  definition  of  the  Hankel  transform,  Zemanian  (1968)  proves  some 
useful  properties  in  his  Theorem  5.5-2.  Although  his  results  are  applicable 
for  any  ’  they  specialized  here  to  sine  and  cosine  functions, 

Jji  =  and  /A=  -  respectively. 

By  Theorem  5. 5-2(5)  with  p  =— ~  ,  one  finds 


With  p  =  +  *2  ,  one  obtains 


J  TOO 

•  *3y  J  dx 

■  -  -j  ,  one  finds 

(C  V 

-  stFtty)  dx 

/*  =  ~  J  ,  one  finds 

/  r<» 

7j7j 

(c.  iv 
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In  each  case,  if  the  integral  on  the  left  does  not  exist  in  the  conventional 
sense,  a  generalized  function  can  often  be  defined  for  it  through  the 
expression  on  the  right. 

For  example,  consider  the  integral 


foffflx)  tor  (Lx)  fa 

X 


Standard  integral  tables,  such  as  Dwight  (1961),  give  the  value  of  the  integral 
to  be  infinity.  If  we  take  the  derivative  with  respect  to  b  and  use  (6.8),  we 
have  0Q 

°  Jo 

\ 

i  r03 

=•  “  z  J  4oc  c/'x  (£|7j) 

As  Butkov  (1968,  p.  312)  points  out,  these  integrals  are  divergent.  If  they 
are  considered  to  be  generalized  functions,  however,  Butkov  states  without 


00 

~  —  \  C«»S  H'lC  Sin  k*  clx 

'O 

OD 


proof  that 


\  kv-  - 

yJn  IN 


as  an 


V 

This  can  be  proven  if  V  =  vo  is  interpretted 

inverse  Fourier  transform.  Thus,  our  integral  becomes 

1  - - 1  (UV 


_1 


2-  Tk 


a^b7- 


Thus,  since 


where  C  is  an  arbitrary  constant.  This  agrees  with  the  generalized  value  for 
this  integral  given  by  Khadem  and  Keer  (1974) . 

Another  integral  which  we  will  have  to  evaluate  is 

'OO  Co 


V 

-/n 


SIH  tx  Sill 


UK  dy  =  (% 


sin  4(x 


J, 


OHM 


tx  d)c 


(uo 


a 


Using  (6.8),  we  find 

OD 


Sin  tx 


sin  ax 


d*  -  — 


At 


ysCO 

ax  be  4 

Jfl 


=  -TTrl+^-^  >o  ($■!?) 


Thus,  we  have 


pC X) 

j  fiv'  tit  tin  UK  d)t  =  E-  %(U~t) 


(6.2°) 


which  agrees  with  Stakgold  (1967,  vol.  I,  p.  286)  and  Davis  (1974).  In  the 
next  section,  the  dual  integral  equation  approach  will  be  used  to  solve  the 
problem  for  the  perfectly  conducting  case  Zg  =  0. 


VII.  QUASI-STATIC  DUAL  INTEGRAL  EQUATION  SOLUTION  FOR  A  SLIT  IN  A  PERFECTLY 
CONDUCTING  PLANE 

The  problem  of  a  slit  in  a  perfectly  conducting  plane  is  of  interest 
because  for  a  typical  graphite  composite  sheet,  Zg  =  0.045  ohms.  This 
means  that  one  can  consider  Zg  =  0  for  some  applications.  Upon  setting 
Zg  *  0  in  (5.31),  one  has  the  dual  integral  equations 
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✓*00 

j  p  (X)  s/m  X? 


'0 

*0D 


f  f(V)  !m  Xf  <i\  =  o 

J6 


OC  <  i> 

•H?  t 


(7.1=) 

0'iy 


A  similar  pair  of  equations  is  solved  in  another  way  by  Clemmow  (1966,  p.  91). 
The  solution  to  these  dual  integral  equations  is  given  by  Sneddon  (1966,  p.  103) 


and  is 


F(X')=  kJ,(H) 


Substituting  this  in  equation  (5.28),  one  obtains 

o“  *  X7  e  ^ 

w.iA  <?+Xx 


H.bfj.CXb) 


-zw  o? 


-X  7  o 


(7.Z) 


Oc  <o 


Jo 


To  evaluate  this  explicitly,  one  needs  to  evaluate  integrals  of  the  type 

SI^X-Z  e  ^  (7,3) 


[ 


CO 

0 


Integrals  of  this  type  have  been  evaluated  by  George  (1962).  The  result  is 

JX 


.  jl  o  l.  l  i  it  o  type  nave  uccu  cvaxuatcu  u  j 

r«  /v  ,  .  -Xwi  i\ 

jG  3,  (Xt)  SIM  X*  e  "X  ~ 


2  *<»*(?)  \ 
=  V  "  V 


(7.4) 


Note  that  for  x=0,  this  becomes 

IV- 2*1  -C^-gTT  (  _ 


V> 


f2' 


59^(2) 


£ 

T 


2C  t> 


(.7.5) 


iiMi 
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If  the  expression  for  z  b  is  multiplied  and  divided  by  (?+  ,  we  find 

the  same  result  as  that  given  by  Magnus  and  Oberhettinger  (1949,  p.  36). 

,c0  JI  f  I 

3,  (11)  />ih  \z  4  c  i  t  (?■  L> 

I  b 


1 -+  J 


The  potential  becomes 


<£  =  -|40  )  1  +•  *)»(**) 


\|0(Ml-2*)‘  +  4x2zz  1 


(7' 7; 


The  x  component  of  the  H-field  is  given  by  equation  (5.9). 

Hx  =  --f£  (7.?) 

After  performing  this  differentiation  and  simplifying,  one  obtains 


=  H0z 


-Ho  ^^bcl 


r 


_  'yJCxV^-Z2/  ^Ayz27-  +  t  ?*) 
2  Oc^b1-?2')2'  4  4  X7-?* 


j  JCxN-  bv-  ?  V*  f  4X7-  i1  ~  *-bu  -  ?r) 
2  "  (Kz  21)*  +  4  xz 


2. 


0-v 


Special  cases  of  interest  are  x=0,  z=0. 


U*  (x,  o)  =  o 


0 


set 

2  ?  t 


(7.109) 


(7. 10  y 
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(7.11) 


The  z-component  of  the  H  field  is  given  by 


2  97 

Performing  this  differentiation,  we  have 


+  9jMWH0|7| 


l  ■JOcHb’-z1)*  +  0£+£-J^) 


1  (x1  +  i>1-  H1)2-  +  4xlZl 

(l  Jo^b1--  rf  rO 


X 

z 


L 


A_ 

■x. 


For  the  special  cases  x=0  and  z=0,  one  obtains 

«•('  +  w^) 


\\i(o,i)=  H 


("7.  n) 


N  <  [>  (J.  I4«) 


VA  o  (tc-  o~  ?)  - 


"  ^  +  {&) 


(7.  14 1>7 


\?l>  l 


[zmTO-21 

)?i 


-  W.(l 


,>i  C7-  k  c; 
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The  equivalent  current  in  the  conductor  can  be  written  as 


=  -2  Ho 


21-f- 


(7./S-; 


VIII.  DUAL  INTEGRAL  EQUATION  FORMULATIONS  FOR  A  SLIT  IN  AN  IMPEDANCE  SHEET 
The  fundamental  quasi-static  dual  integral  equation  for  the  problem 
of  a  slit  in  an  impedance  sheet  is  given  by  (5.31).  As  it  stands,  (5.31a) 
is  not  in  a  form  whose  solution  appears  in  Sneddon  (1966).  There  are  several 
dual  integral  equation  pairs  that  can  be  obtained  from  (5.31).  The  first 
results  from  direct  differentiation  in  (5.31a).  Using  (6.10),  we  have 

SfeXz  ^  l*Kl,  {%\e) 

-225X  +  sin  Xi  dx  —  °  (?-IW 

If  the  methods  of  Sneddon  (1966)  are  used  to  reduce  this  to  a  different  form 
for  solution,  one  finds  that  the  normalized  integral  equation  to  be  solved  Is 


of  the  form 


ivkw  *  -J j  l.to  « i  -  i }j=7==r-j  J« 

f1  .  .  C  y(ZfcF)' 

-  _  . . ,  is_  ill  a*  i-  u  4-  1  M  iu 

(if  J  du  du  ^  (uT fuTT  , 


K  is  the  elliptic  integral.  In  order  to  derive  this  result,  heavy  use  of  the 
generalized  Hankel  transformation  had  to  be  made.  Most  integrals  that  arise 
have  value  infinity  when  treated  in  the  conventional  sense.  This  equation 
appears  to  be  too  complicated  to  solve  easily,  and  therefore,  the  dual  integral 
above  will  be  rejected. 

Another  dual  integral  equation  can  be  obtained  from  (5.31)  by 
manipulating  the  second  integral  in  (5.31a).  From  (6.10),  one  has 


I*1  ( 60 

^  K^-Tcv)  sii*(xy)  <!x  ~  ~ 


By  identifying  a  new  function  F(x)  ■  x  f(x),  one  finds  that 

CfCti  si*  \i  d\  -  -  -j?  1 F  sIm  x  2 

</*  'o 


Substituting  this  for  the  second  integral  in  (5.31a),  one  obtains 


(s- 4; 


*"  Xl  T"  Go  t2 2s 

/» 

l  F(W  ^ 

/•l 


(?.5) 


JX  l  - 


-  o 


22s  -  si*\* 


=  (P£ fV  +  As^B 


+  C  s<jwfe) 


1*1  >b 

(?y 


where  A,  B,  and  C  are  constants  of  integration.  The  C  sgn(z)  term  is 
included  because  the  limits  l3(>  \>  only  allow  such  a  term.  Indeed,  any  term 
which  has  C  sgn(z)  behavior  for  and  any  behavior  whatsoever  for  |?l<t  is 


allowed.  The  dual  integral  equation  becomes 


[  m  *•>.  ^  ll"1 

fftsfa  -  i  =  [a  *  (%%£)}  -18 


?„+22. 


la  I  <  b  (Z^) 


(?.%) 


-vC^fe)  lw>t 


This  equation  will  be  converted  to  a  Fredholm  equation  of  the  second  kind  in  the 
next  section. 


The  constants  A,  B,  and  C  in  equation  (8.7b)  need  to  be  evaluated. 

The  evaluation  of  these  constants  is  not  straightforward  and  a  unique 
result  is  not  forthcoming.  Values  for  A,  B,  and  C  will  be  assumed  here 
based  on  certain  assumptions.  Since  this  problem  has  apparently  not  been 
solved  before,  there  are  no  results  to  compare  the  results  of  this  work  to. 

This  makes  verification  of  the  assumptions  used  to  find  A,  B,  and  C  impossible 
at  this  time. 

To  find  a  value  for  B,  one  can  assume  that  since  the  left-hand  side 
of  equation  (8.7b)  is  an  odd  function  of  z,  the  right-hand  side  must  also  be 
an  odd  function  of  z.  Therefore,  the  constant  B  must  be  zero 

The  assumption  that  will  be  used  to  find  a  value  for  A  is  that  the 
right-hand  side  should  remain  finite  as  z  approaches  infinity,  thereby  making 

2 

4  ■  - 2  zs  <  xA r  >  <8-8> 

0  s 

Note  that  his  is  not  the  only  assumption  that  could  be  made.  If  A,  B,  and  C 
were  all  taken  to  be  zero,  for  instance,  the  right-hand  side  would  be 
proportional  to  z.  If  equation  (8.7)  were  considered  to  be  a  conventional 
dual  integral  equation,  this  would  almost  certainly  be  disallowed  because 
the  right-hand  side  would  blow  up  at  infinity.  If,  on  the  other  hand,  we 
recall  that  we  are  allowing  generalized  function  solutions,  then  since  z 
is  a  distribution  of  slow  growth  over  the  interval  (1,  00  ) ,  a  right-hand 
side  proportional  to  z  is  not  disallowed  in  the  generalized  sense. 

There  are  at  least  two  reasons  to  assume  that  A  is  given  by  (8.8). 

First,  we  formulated  the  problem  by  using  an  unknown  function  which  is 
assumed  to  be  small  far  from  the  slit.  The  fields  far  from  the  slit  were 
already  accounted  for.  Second,  it  gives  agreement  for  the  perfectly  conducting 
case.  Let  Zg=0  in  equation  (8.7b).  Then  we  obtain 

^  =  A*  121  ?l 

where  B  has  been  taken  to  be  zero.  For  the  perfectly  conducting  case,  we  found 
(Section  VII)  that  F(  X  )  =  b  Jj(Xb).  Substituting  this  in  the  above  integral, 
one  obtains 


Thus,  one  has 


-AIa)JS 


A? 


i^i 


If  we  let  z  get  large,  we  see  that  A=0.  This  agrees  with  our  assumption 
that  A  =  -2ZsZq/(Zq+2Zs)  because  Zg=0  In  this  case.  Assuming  that  (8.8)  for 
A  is  correct,  the  dual  integral  equation  (8.7)  becomes 


c  1*2/ 


C  will  be  evaluated  later.  This  equation  will  be  used  later. 

There  is  another  dual  integral  equation  that  one  can  obtain  from 
(5.31)  by  using  (6.8)  through  (6.11).  Differentiating  (5.31b)  and  using 
(6.11)  with  f(x)  =  F(x)/x,  one  finds 

f°°  ~Z 

\  F()0  cos(X?)<iX  -  fetch 

'Jo  0 


By  using  (6.8)  and  (6.11),  (5.31a)  can  be  written  as 

JT0)[22S  -  +  C0  l 

where  Cq  is  an  arbitrary  constant  of  integration.  If  we  assume  that  Cq  is 
not  proportional  to  Zg,  then  for  the  perfectly  conducting  case  when  Zg=0, 
F(  X  )  =  b  Jx(b  X)  and 

-  C 

Jo 


Since 
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(  ‘osCW-^ 

-'a 


A  _ 


i-G)' 

o 


IVK  l 


111  >  \> 


one  has  Cg=0.  The  dual  integral  equations  for  this  case  then  become 

r00 

\  f(x>  £os  x?  A  =  jTTzzs  |?| 

^f(X)[223  -  X-z  A  =  Z4  +  2 ?s 

A 

If  the  techniques  of  Sneddon  (1966)  are  applied  to  these  equations,  one  finds 
that  an  auxiliary  integral  equation  that  needs  to  be  solved  is 

22,  Wt)  +  j 

r00 


2b  / 


L 

^  "^TT  V 


(■ 


Zb 

T  V?0  +  22s 


i?/  7>t> 


+c 


where  Cj  and  are  constants  of  integration.  In  order  to  derive  this  result, 
considerable  use  of  the  generalized  Hankel  transformation  of  Zemanian  (1968) 
had  to  be  made.  The  kernel  turns  out  to  be  (6.17)  which  is  infinity  if 
interpreted  in  the  conventional  sense.  Grosjean  (1972)  briefly  presents 
some  identities  for  the  $n  |t^  -  z^\  kernel.  Because  of  the  fact  that 
conventional  mathematics  breaks  down  in  this  case,  the  above  dual  integral 
equation  will  be  avoided  for  further  study.  The  dual  integral  equation 
(8.12)  is  the  only  equation  we  will  choose  for  further  study.  The  next 


section  presents  a  method  for  changing  this  dual  integral  equation  into  a 
Fredholm  integral  equation  of  the  second  kind. 


A.  FORMULATION  FOR  LARGER  VALUES  OF  SHEET  IMPEDANCE 
The  dual  integral  equation  of  (8.13) 

(  F(»  4-  ~  *  tel  <t 

Jo  A 

^  C  ft) ^l>t  (?*IW 

give  the  best  results  when  Z  is  large.  C  will  be  assumed  to  be  zero  for 

s 


convenience.  For  infinitely 

large  Z  ,  we  have 

C  f(y)  si*  x^ 

4-  =  ° 

zcb 

(?.2«) 

“'O 

(€0 

)  F(X)  X? 

Jo 

4=0 

(p.zy) 

or  F(  \  )  =  0.  This  makes 

*4 

04  \?o 

0  "  %  CO 

(?v 

equal  to  -H^z,  all  x,  as  we  expect.  The  dual  integral  equation  ,  in  this  form, 
gives  us  a  result  which  has  F(^)  smaller  and  smaller  for  larger  and  larger  sheet 
impedances. 

The  dual  Integral  equation  (9.1)  can  be  reduced  to  a  single  Fredholm 
integral  equation  of  the  second  kind  in  a  straight-forward  manner  if  we  let 
C=0.  Let 
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(9.4) 


F(^)  =  G(\)  +  H(\) 

where  G{\)  is  a  function  which  satisfies  (9.1a)  and  H(  X)  is  such  that  the 
integral  in  (9.1a)  is  equal  to  zero,  that  is, 

^  sfn  Xs?  ~~  ®  ^ i=»  (*  ^ 

Sneddoti  (1966,  p.113)  gives  formulas  for  determining  G(^)  and  H(/\  )  •  In  the 
present  case,  his  formulas  give 

<?■« 

,  60 

( 

where  h(t)  is  an  unknown  function  which  will  be  solved  for  later.  By  sub¬ 
stituting  (9.4)  into  the  dual  integral  equation  (9.1),  one  obtains 


The  integrals  involving  G(  X)  on  the  right-hand  sides  can  be  found  in  standard 
integral  tables.  They  are 
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(?'\0) 


where  is  the  Legendre  function  of  the  second  kind  as  defined  by  Magnus, 

Oberhettinger,  and  Soni  (1966,  p.  176).  In  terms  of  more  familiar  functions. 


&(»>  -  =  -fj*  +  hr  &  (|rf)]  2  7 1 


Substituting  these  integrals  of  G(X)  in  (9.8),  and  simplifying  yields  the  new 
dual  integral  equation 
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^H(X)  A?  y-  =  O  (9.|2sj 


*H  (X)  [ 


29c  - 


’]  «- 


U1  >  L 


Substituting  the  expression  {9-0  for  H(X)  in  the  first  equation,  one  finds 

j  HQo  fl-Xs  4  = 

^  j_  r°0  ,oo 

-(?T^)  \  Sr^(L^t)  Uu><(d  siu  X?  d\  cR 

x  ,r'  ,<* 

=  (~f 'C  \  Mo  smxtj  sih  X*  JX  J-t  (?■  n) 
J) 

This  assumes  that  h(t)  is  such  that  the  interchange  of  the  order  of  integration 
is  permissible.  By  using  the  theory  of  distributions,  one  can  show  that 

si*  X**  dV  —  "z~  S  *™  ”2; )  (9  I4J 

-'A 


Substituting  this  in  the  above  integral,  one  obtains 

$TlOO  r*X,  ^  - 


3S 


o 


1  .  J. 


(9.10 


(z~)7’  t  1  iidj 


This  shows  that  such  an  H(,A*)  satisfies  (9.12a). 

Substituting  (9.7)  in  the  second  equation,  one  has 


(?■  IQ 


Again  assuming  that  h(t)  is  such  that  the  order  of  integration  can  be  changed, 
one  has 


j  [2ZS  -  s'-4fA) 


S(h  ? 


\  Jt  = 


Zo  \fiY-  l' 


but 


^  sikUX  *  X  JX  - 

-Z^sX^-2)  -  iJL 


(?.n) 


(?.  is) 
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Substituting  this  in  the  above  equation,  one  obtains 


Using  the  sifting  property  of  the  delta  function  in  (9.15),  one  obtains  the 
Fredholm  integral  equation  of  the  second  kind 


Using  the  change  of  variables  /T'  =  bt  in  the  integral  and  GJl*  =  kZ^,  °ne  finds 


By  using  the  fact  that 
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the  integral  equation  can  also  be  written  as 

r  od 


(?.  zz) 


The  dual  integral  equation  has  now  been  reduced  to  a  single  Fredholm  integral 

equation  of  the  second  kind.  If  this  equation  can  be  solved  for  h(  —  ),  then 

b 

the  solution  to  the  dual  equation  (9.1)  can  be  written  as 


In  order  to  show  the  necessity  of  using  the  theory  of  distributions,  it 
is  helpful  to  point  out  what  happens  if  (8.4)  is  not  used  to  simplify  (5.31). 
Equation  (5.31a)  can  be  written  as 


22, 


/><*  rOO 

i_—\  F(X)  ~~y  +- wy \  ^  ~~  o  (9,25) 

*  4)  Jo 


Substituting  (9.4)  in  the  above,  one  finds 
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00 


Jo 


(r(\)  SUo X?  dx  (fzQ 
z?L 


From  (9.9),  the  first  integral  on  the  right-hand  side  is  equal  to  zero.  The 
second  integral  can  be  reduced  to 


'CO 


lUUf  dX  c^v 


This  integral  can  be  expressed  in  the  form  of  a  Weber-Schafheitlin  integral  as 
given  by  Watson  (1944,  p.  398-404).  Unfortunately,  using  Watson's  criterion, 
the  convergence  of  this  integral  is  not  assured.  Indeed,  it  is  not  listed  in 
any  standard  integral  tables.  It  would  thus  appear  that  the  dual  integral 
equation  approach  can  not  be  used  to  solve  this  problem  using  classical  math¬ 
ematics.  We  will  assume  that  all  of  the  steps  leading  to  the  integral  equation 
(9.23)  are  justifiable.  The  justification  rests  on  a  foundation  of  the 
theory  of  distributions. 

B.  FORMULATION  FOR  SMALLER  VALUES  OF  SHEET  IMPEDANCE 

Since  for  the  composite  problem  at  hand,  we  know  that  Zg  is  very  small, 
typically  0.045  ohms  and  that  the  above  formulation  is  best  for  large  Zg,  it 
is  necessary  to  find  a  formulation  good  for  small  values  of  sheet  impedance. 

One  possible  way  to  do  this  would  be  to  subtract  the  known  solution  bJ^(\b) 
for  the  perfectly  conducting  case  Zg  =  0  from  the  unknown  function  F(X)  of 
(9.1)  and  try  formulating  a  dual  integral  equation  in  terms  of 


Fg(  X)  =  F(X)  -  bJ^bA) 


(9.28) 
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where  the  s  subscript  means  small  values  of  Z  .  Adding  and  subtracting  terms 
of  bJ^(bX)  to  F(\)  in  (9.1)  gives  the  dual  integral  equations 


^Fs0) At  ^  = 


f 


00 


Zo+22* 


c«*o 


5oFs(X){22s-  si«  X?  A.  = 

00  , 


v  C 


Erdelyi,  et  al  (1954,  p.^9)  have  evaluated  the  integrals  on  the  right-hand  sides. 

2  ^L> 


J|  0*\)  ^  "  J  ^ 


121  -t’/iM?’ 


2> 


+T'^”'(t)  icl 


b\3|CWciVi\?  — 


TT 


Substituting  these  results  in  (9.29)  gives  the  dual  integral  equations 
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^  F^()0  _  —  — 

Jq  *  2 .Q 


2?< 


+  2  2, 


?<l  (?.  32a) 


-  _  ALJ/J  7  •  l  iX_ 

2zs~  —js,,a  x*  x 


*420irtL 


4 


*  C 


z?s  w+fcrrjp 
2>t 


Proceeding  In  the  same  manner  as  before,  one  finds  that 


B(»=  G,(X)HlsO0  =  -(ir 


"  4  (U)  +■ 


+ 


L-AJ  *in(Li:  X)lof  (VJt  &?}) 


/jX  225  ^ 

^2'  x4?o  i?i  +  jzv-  kT 


(9.34; 
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If  C  is  chosen  to  be 


C  =  -  ikZQ^|-  (9.35) 


then  one  has  to  solve  the  integral  equation 


It  Is  not  known  if  the  above  value  for  C  is  correct.  It  was  chosen  to  allow 
an  easy  solution  in  terms  of  Legendre  polynomials. 


X.  EXPRESSIONS  FOR  POTENTIALS  AND  FIELDS 

The  unknown  function  F(A)  in  the  dual  integral  equation  (8.12)  has  in 
the  last  section  been  specified  in  terms  of  a  known  function  G(\)  and  an 
unknown  function  H(  X  )  .  It  is  of  interest  to  find  the  potentials  and  fields 
in  all  space  once  F(X)  has  been  found.  For  the  formulation  of  Section  IX. A., 
the  potential  is  given  by  equation  (5.28)  with  F(X)  replaced  by  G(A)  +  H(X)* 


-H. [[600+ 


(lO.laJ 
(10. 1 l) 
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G(  X)  and  H(X)  are  given  by  (9.6)  and  (9.7),  respectively.  For  the  formulation 
of  Section  IX. B.,  the  potentials  are  again  given  by  (5.28)  with  F(  X)  given 
by  equation  (9.28) 


F(y>=  fjx)  -t- 1  T,(U)  (io.2) 


where  F  (X)  =  G  (X)  +  H  (X)  is  given  by  (9.33).  In  both  formulations, 
s  s  s 

one  has  to  do  potential  integrals  involving  functions  of  the  form  G^(X)  = 
b(bX  )^3/2^  X )  •  Substituting  this  in  (10. 1),  one  finds  integrals  of  the 
form  ^Gradshteyn  and  Ryzhik  (1980,  pp.  491-492)J 


4  -  fx> 


-Xlxi 


A 

X 


=  itCuh^U)  AAi  e  ' Xlx'-^  = 

Jo 


0°*) 


where  if  ^  (z^  -  b^)  then  s=l  else  s=0.  The  H(/\)  and  Hs(\)  contribut 
can  be  reduced  to  integrals  of  the  form  f^  and  f^  where 

•CO 


ions 


^  =  H(X)  si*Xe  e  X  ^  ^ 


^  ^ j1"  1  \  ^  siv\(bX£)k(W  dt 


z  e 


-X,y|  JA  _ 
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The  inside  integral  can  be  evaluated  analytically  |j:rdelyi  et  al  (1954,  p.  159)J 


si*\A3-  e  d\  — 


_ 2ixi  ?yt 

t*1  -*•  (bt  -  +  (u+?fj 


C/O.s) 


If  for  x=0,  this  integral  is  interpretted  as  a  distribution,  then 

/~co  —  Alvl 

\  si* <?  d\  —  for  X=  O  (/°£>) 

Jo 

where  ^  is  the  delta  distribution.  The  integral  f  then  becomes 

j_  rco 

fH  is  given  by  this  integral  with  h(t)  replaced  by  h  (t).  For  the  special 


case  x=0,  one  obtains 

~  _ 4-  .J-f00 


It  is  interesting  to  note  that  h(t)  is  just  the  value  of  f  evaluated  on  the 

n 

impedance  sheet.  The  final  potential  integral  that  needs  to  be  evaluated  is 


This  integral  has  already  been  evaluated  in  Section  VII. 
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The  complete  expressions  for  the  potential  when  using  the  formulation 
of  Section  IX. A.  is 


The  potential  when  using  the  formulation  of  Section  IX. B.  is 


(Vo.  I  ok) 


where  PEC  is  the  abbreviation  for  Perfect  Electrical  Conductor. 
The  x  and  z  components  of  the  H-field  may  be  obtained  from 


Go.  it) 


XI.  REDUCTION  OF  THE  FREDHOLM  INTEGRAL  EQUATIONS  TO  INFINITE  SYSTEMS  OF 
LINEAR  ALGEBRAIC  EQUATIONS  IN  AN  INFINITE  NUMBER  OF  UNKNOWNS 
Fredholm  integral  equations  similar  to  (9.20)  and  (9.36)  have  often  been 
solved  by  expanding  the  unknown  and  the  right-hand  side  in  series  of 
orthogonal  polynomials  such  as  Chebyshev  or  Legendre  polynomials.  Morar’ 
and  Popov  (1970),  Popov  (1972),  and  earlier  untranslated  works  of  Popov  (in 
Russian)  use  this  solution  procedure.  England  and  Shall  (1977)  use  the 
technique  for  solving  some  elasticity  problems.  Aleksandrov  and  Kovalenko 
(1977)  also  solve  elasticity  problems. 
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A.  FORMULATION  FOR  LARGER  VALUES  OF  SHEET  IMPEDANCE 

The  Fredholm  integral  equation  that  needs  to  be  solved  for  this  case  is 
given  by  (9.20).  Substituting  z=ub,  one  obtains 

{Tj-  - 

By  using  symmetry  (the  right-hand  side  is  an  odd  function)  and  two  changes  of 
variables,  one  obtains  the  equation 


O',  i) 


-b 


( IU ) 


*  —2  1 

Let  h  (t)  *»  t  h(  —  ).  Then  the  above  equation  becomes 

r  I 


2TT^5  r-1. 

- IS  T  —  -  3 


^cs)  4-  r  t,%)A(s-vi  jv = 

n/ 


0  !■!) 


-b  l?  /  v 


?o  \± 
2o  +22j  j  S’- 


i  S-f 


The  right-hand  side  can  be  expressed  as  a  power  series. 


r 


L 


2  H  4-  I 

+  i-t ( }  z £o (zn+ dC2m+ 


3) 


Cii.  4; 


This  is  an  odd  function  and  can  also  be  expanded  in  odd  order  Chebyshev 
polynomials  and  Legendre  polynomials  j.  Only  the  Legendre  poly¬ 

nomial  solution  will  be  attempted  here  because  the  Chebyshev  solution 
apparently  gives  fields  which  are  more  singular  than  they  should  be  near  the 
edge.  Since  by  Luke  (1969,  vol.  I,  p.  277) 


f"  --  2  Oi'^J 

Y*  -  0 

o  ,  2^+3. 

(h-w)! 


The  power  series  can  be  written  in  terms  of  Legendre  polynomials. 

/2-H  +  l  CO  , 

°°  t  «=c  I 

22 


£0  (2*+>)0+3) 


It  can  easily  be  shown  that 

CTO  Y\ 


2 


*  \ 

£t0  (zn+oo+i)  irL  W 

(ii. t) 


h-  o 


lv»t=o 


(n: 1) 


Using  this  identity,  one  finds 


Caw 
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where 


OO 


2*f  T- 


1  -  T  CZmHJltWm-1)!  A |  gy 

^  I  (t^+0(^lM+3)[2(’^fM  +l)]j  (fc-n  +l)J 

* 

Assume  that  the  unknown  h  (t)  can  be  expanded  in  a  Legendre  polynomial 
series. 


h*(t) 


CO 

■  2  bnp2n_,(t) 

n=l 


(11.9) 


2  * 


In  order  to  use  orthogonality,  t  h  (t)  must  be  expressed  as  a  series  of  Legendre 
polynomials.  Identity  8.915  5.,  page  1026  of  Gradshteyn  and  Ryzhik  (1980)  can 
be  used  to  do  this.  One  finds  that 

CH./o) 


a0U)  = 


(zM(zt) 


a 


(<?l hlW-l) 

Zi  (2{-l) 


I  +  2 


I) 


n  rlt)  -  H-W-V 

iW~  cu-m-3) 


(II.  Ill) 


(ii.IIW 


(II.  He) 


Note  that  a^fl) 


0.  Substituting  (11.8),  (11.9),  and  (11.10)  in  (11.3),  one 


z _ 

0*  y  \ t-n- -  hjl  -  f| 


(n.&) 


Multiplying  through  by  P2m-l^s^’  inte8ratin§  from  -1  to  1,  and  simplifying, 
one  obtains 

<30  _  _  Z 


*• ;) J 


=  I11  (ir  (■ 


^o+2?s. 


VHff-l 

(.IU4) 


This  is  an  infinite  set  of  equations  in  an  infinite  number  of  unknowns  which 
must  be  truncated  and  solved  on  the  computer.  This  will  be  done  in  a  later 
section. 

B .  FORMULATION  FOR  SMALLER  VALUES  OF  SHEET  IMPEDANCE 

The  Fredholm  integral  equation  that  needs  to  be  solved  for  this  case  is 
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given  by  (9.36).  Substituting  z=ub  and  T  =tb  in  this  equation  gives  the 
equation 


=  -t 


\(-L\X 


(11.15) 


I  ^(SLY-  vp(M> 

■ill'Zo  lUl+Jui-i 


|m|>| 


The  last  term  in  this  equation  can  be  simplified  if  numerator  and  denominator 
multiplied  by  \U\  —  ,  •  ^ne 


are 


I u\  ±[U^T 

By  using  symmetry  and  two  changes  of  variables,  one  obtains 


(ii.it.) 


=-^f{ 


Vz.[7-\tl  _\JL  \  %  -f 

?o+2?sJl 
i- 


^KUrl] 


^n> 

Let  h*(t)  =  t“2h  (  -  ).  Substituting  this  in  the  above  equation,  one  obtains 
s  s  t 
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2tt2s  ,t-I,V)  , 

TIEUT  ^(s;  + 


f 


dv/  = 


o-w  «■ 

Except  for  the  last  term  on  the  right-hand  side,  this  equation  is  similar  to 
(11.3).  By  using  Luke  (1969,  vol.  I,  p.  38),  the  last  term  becomes 

-  til 


(it.  i?; 


Itl<  I 


(  )  is  the  Pochhammer's  Symbol.  Substituting  (11.5)  in  this  equation  and 

using  (11.7),  one  obtains 

I- d-t'^5  ={f  e,?26_,(H  Quo.) 

{-( 


where 


eA  = 


C k+4-bi)!  (z{-V  z 
(in-i7Ty 


:f*\ 
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Expanding  hg(t)  in  a  Legendre  polynomial  series  and  proceeding  as  before,  one 
obtains  the  infinite  set  of  equations  in  an  infinite  number  of  unknowns 

pr-1)  i,  J-A  = 

✓  7 

Note  that  if  Zg=0,  then  bn=0  for  n=l,2,...^®  as  expected. 


I ) 

■Ob 
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XII.  NUMERICAL  SOLUTION  OF  EQUATIONS  (11.14)  and  (11.21) 

The  simultaneous  equations  (11.14)  and  (11.21)  were  solved  on  the  computer 
by  truncating  the  limit  on  the  summations  to  N,  a  variable  in  the  computer 
programs.  Doing  this  for  (11.14)  yields  the  result 


presents  a  computational  difficulty.  The  infinite  series  (11.8b)  turns 
out  to  be  very  slowly  converging.  Let 

J~  CD 


where 


= 


C^vvi+I)[ 


(2w  +?)  [2(^+0 +i;J/  (i*i-n +o! 


(12.1; 


The  next  term  in  the  series  can  be  computed  from 


The  first  term  in  the  d  series  is  given  by  D  , (n) 

n  b  J  n- 1 

next  series  (d  ,.)  is  given  by  D  (n+1). 

n+i  n 


Dh0"+|) 


(Zh-l)  h  (2^-H  I) 


.  The  first  term  in 


the 


(12.5V 
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The  numerical  values  for  the  dn's  are  given  in  Table  12.1  along  with  the  upper 
index  of  summation  I  in  (12.2).  As  can  be  seen,  a  large  number  of  terms 


Table  12.1.  Computed  values  of 


n 

I 

d 

n 

1 

209 

0.7988899796 

2 

523 

0. 1205192866 

3 

892 

0.0382577779 

4 

1308 

0.0165125030 

5 

1765 

0.0085110354 

6 

2256 

0.0049279487 

7 

2776 

0.0030974268 

8 

3324 

0.0020692245 

9 

3896 

0.0014488469 

10 

4491 

0.0010530191 

11 

5106 

0.0007888795 

12 

5740 

0.0006060148 

13 

6393 

0.0004754722 

14 

7062 

0.0003798273 

15 

7748 

0.0003081716 

16 

8449 

0.0002534407 

17 

9165 

0.0002109242 

18 

9895 

0.0001773995 

19 

10638 

0.0001506120 

20 

11395 

0.0001289524 

21 

12164 

0.0001 1 12509 

22 

12945 

0.0000966439 

23 

13737 

0.0000844842 

24 

14000 

0.0000742759 

25 

14000 

0.0000656441 

26 

14000 

0. 0000582972 

27 

14000 

0.0000520049 

28 

14000 

0.0000465846 

are  required  to  give  seven  digit  accuracy.  If  more  than  14000  terms  are  needed, 
the  summation  is  truncated  there.  As  a  check,  the  value  of  the  Legendre  series 
in  (11.8a)  was  evaluated  at  ten  equally  spaced  points  between  zero  and  one  and 
compared  to  the  values  obtained  by  other  methods.  The  two  values  so  obtained 
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matched  for  all  points  to  approximately  four  digits  of  accuracy.  While  this 

was  only  a  spot  check,  it  did  show  that  the  Legendre  series  in  (11.8a)  is  a 

proper  representation  for  the  function  on  the  left  in  (11.8a),  and  that  the 

d  's  in  Table  12.1  are  correct  to  at  least  four  decimal  places, 
n 

The  simultaneous  equations  (12.1)  were  solved  on  the  computer  for  various 

matrix  dimensions  N,  wavelengths  X,  slit  half-widths  b,  and  sheet  impedances 

Z  .  It  was  hoped  that  b,  would  be  much  larger  than  the  other  b  ’  s  and  so 
si  n 

would  provide  a  one  term  approximation.  Unfortunately,  this  was  not  the  case. 

For  (A,b,Z  )  =  (3m, 0.5cm,  0.045)  and  N  from  10  to  28,  most  of  the  b  's  were  of 

S  -4-6  n 

the  same  order  of  magnitude  (10  to  10  ).  Convergence  was  quite  good  in  that 

the  values  of  the  b^'s  did  not  vary  much  with  changing  N.  For  this  particular 

choice  of  (A,b,Zg),  it  does  not  appear  that  a  one  term  solution  could  be  used 

to  give  suitable  results.  One  thing  should  be  kept  in  mind  though.  This  is 

simply  the  possibility  that  in  equation  (10.9),  f„  could  be  significantly 

H 

lower  in  magnitude  than  the  f_  term.  As  they  stand,  the  equations  are  in 
such  a  form  that  it  is  difficult  to  evaluate  the  relative  importance  of  the 
fG  and  f^  terms.  The  individual  terms  when  h(t)  in  (10.7)  is  replaced  by  the 
Legendre  series  are  of  the  form 


A  general  closed  form  expression  for  any  n  for  this  integral  has  not  yet  been 
found. 

The  wavelength  X  was  made  larger  to  X  =  30m.  The  solution  for  this  case, 

(X,  b,  Zg)  =  (30m,  0.5cm,  0.045),  is  given  in  Table  12.2.  From  the  results, 

one  can  see  that  only  the  first  and  second  terras  probably  need  to  be  retained. 

The  cases  (30m,  0.05cm,  0.045)  and  (3m,  0.5cm,  0.45)  were  similar  in  that 

the  first  three  terms  were  dominant.  These  three  cases  were  fairly  special, 

however.  In  cases  such  as  (3m,  0.5cm,  1,000,000),  (30m,  0.05cm,  0.045), 

(300m,  0.5cm,  0.045),  and  (3m,  0.5cm,  0)  at  least  ten  of  the  b  's  were  of 

-4 

comparable  magnitudes.  The  decimal  values  ranged  from  2.5x10  for  (3m, 

— 18 

0.5cm,  0)  to  4x10  for  (3m,  0.5cm,  1,000,000).  Due  to  the  very  small  size 
of  the  coefficients  for  large  Z  ,  it  could  be  argued  that  f^  could  be  neglected 
in  equation  (10.9). 
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Table  12.2.  Values  of  b 

n 


N=  28  LAMBDA:  30.00000000  B=  0.00500000  ZS= 

0.04500000 


I. 

SOLUTI ON ( I ) 

1 

-0.0001286771 

0.0000523017 

2 

-0.0000292388 

0.0001135824 

3 

0.0000412771 

0.0000374772 

4 

0.0000097743 

-0.0000068537 

5 

-0.0000035659 

0.0000054503 

6 

0.0000039863 

0.0000043682 

7 

0.0000009835 

-0.0000000319 

8 

0.0000000985 

0.0000016002 

9 

0.0000008883 

0.0000007096 

10 

0.0000001430 

0.0000003934 

11 

0.0000002808 

0.0000005342 

12 

0.0000002166 

0.0000002411 

13 

0.0000001093 

0.0000002722 

14 

0.0000001443 

0.0000001977 

15 

0.0000000804 

0.0000001525 

16 

0.0000000823 

0.0000001421 

17 

0.0000000645 

0.0000001045 

18 

0.0000000520 

0.0000000969 

19 

0.0000000478 

0.0000000784 

20 

0.0000000375 

0.0000000678 

21 

0.0000000348 

0.0000000600 

22 

0.0000000280 

0.0000000498 

23 

0.0000000272 

0.0000000461 

24 

0.0000000195 

0.0000000391 

25 

0.0000000237 

0.0000000332 

26 

0.0000000129 

0.0000000363 

27 

0.0000000188 

0.0000000172 

28 

0.0000000156 

0.0000000417 

Equation  (11.21)  has  also  been  solved  on  the  computer.  Truncating  the 
limit  on  the  summation  to  N,  one  obtains 
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The  d^'s  have  been  calculated  and  are  presented  in  Table  12.1.  Only  the  e^'s 
need  to  be  calculated.  Let 


1=00 

e*  =  feMh  Z  Uil 


where 


fa -9) 


r  n\  =  (n+l+O! 

(^<(zM4-0)!(*-4+i)! 


The  next  term  in  the  series  can  be  computed  from 

r  fei* +!)(?-*< +3)  r  „  r 

4(i«  +2) (zu*  +2h IX* - 1  +y  ^  ^  ^ 

The  first  term  in  the  series  is  given  by  E^  ^(k).  The  first  term  in  the 
next  series  (ek+j)  is  given  by  Ek(k+1) .  The  relation  between  the  two  is 


r  ft  n  _  JiiiXiM L_p  (4) 

|)(4^30(44*-l) 


Oz.iuJ 


Eo(l)  _ 


(iz.nk) 


The  numerical  values  for  the  e^'s  are  given  in  Table  12.3  along  with  the  value 
of  the  upper  index  I  in  equation  (12.8).  As  can  be  seen,  the  e^'s  are 
extremely  slowly  convergent. 

Equation  (12.7)  was  solved  on  the  computer  for  various  values  of  matrix 
dimension  N,  wavelength  A,  slit  half-width  b,  and  sheet  impedance  Zg.  It  was 
hoped  that  a  one  term  result  would  be  obtained.  For  (A,b,Zg)=(3m,  0.5cm,  0.045), 
the  first  ten  bn  coefficients  were  of  comparable  amplitude  as  shown  in  Table  12.4. 
The  trend  in  the  data  for  (30m,  0.5cm,  0.045)  is  similar  to  that  in  Table  12.2. 


Table  12.3.  Computed  values  of  e^ 


k 

I 

ek 

1 

532 

0. 6437826847 

2 

1361 

0.1438840258 

3 

2282 

0.0620739102 

il 

3276 

0.0342583375 

5 

4327 

0.0216270308 

6 

5425 

0.0148645553 

7 

6563 

0.0108335385 

8 

7735 

0.0082413892 

9 

8938 

0.0064777034 

10 

10169 

0.0052241416 

1 1 

11425 

0.0043016205 

12 

12705 

0.0036031900 

13 

14000 

0.0030618054 

19 

14000 

0.0026333372 

15 

14000 

0.0022886738 

16 

14000 

0.0020073083 

17 

14000 

0.0017746398 

18 

14000 

0.0015800448 

19 

14000 

0.0014156449 

20 

14000 

0.0012754983 

21 

14000 

0.0011550553 

22 

14000 

0.0010507849 

23 

14000 

0.0009599128 

24 

14000 

0.0008802349 

25 

14000 

0.0008099822 

26 

14000 

0.0007477220 

27 

14000 

0.0006922837 

28 

14000 

0.0006427035 

The  magnitudes  of  the  first  three  coefficients  are  a  factor  of  ten  larger  than 
the  other  coefficients.  For  (3m,  0.5cm,  10  all  the  coefficients  are  of 
the  same  order  of  magnitude  and  were  in  the  neighborhood  of  3x10  It  does 

not  appear  that  it  can  be  predicted  on  the  basis  of  an  analysis  of  the  data 
when  more  or  fewer  coefficients  will  be  needed.  The  relationship  appears 
quite  complex.  It  does  not  appear  that  any  simplification  occurs  when  the 
"small  sheet  impedance"  formulation  of  XI. B.  is  used. 
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Table  12.4.  Values  of  b 

n 

.Ns.  26  LAMBDA:  3.00000000  Bs  0. 00500000  ZS  = 

0.04500000 


I,  SOLUTION ( I ) 


1 

-0.0000004882 

-0.0000052007 

2 

-0.0000022124 

-0.0000072143 

3 

-0.0000040522 

-0.0000067656 

4 

-0.0000052469 

-0.0000051073 

5 

-0.0000055239 

-0.0000030881 

6 

-0.0000050141 

-0.0000013140 

7 

-0.0000040616 

-0.0000000909 

8 

-0.0000030213 

0.0000005533 

9 

-0.0000021297 

0.0000007617 

10 

-0.0000014770 

0.0000007240 

11 

-0.0000010475 

0.0000005925 

12 

-0.0000007802 

0.0000004546 

13 

-0.0000006134 

0.0000003437 

14 

-0.0000005033 

0.0000002635 

15 

-0.0000004246 

0.0000002069 

16 

-0.0000003646 

0.0000001663 

17 

-0.0000003171 

0.0000001360 

18 

-0.0000002784 

0.0000001129 

19 

-0.0000002466 

0.0000000948 

20 

-0.0000002199 

0.0000000804 

21 

-0.0000001974 

0.0000000687 

22 

-0.0000001782 

0.0000000591 

23 

-0.0000001618 

0.0000000512 

24 

-0.0000001477 

0.0000000442 

25 

-0.0000001346 

0.0000000379 

26 

-0.0000001229 

0.0000000354 

27 

-0.0000001246 

0.0000000321 

28 

-0.0000001183  • 

-0.0000000114 

XIII.  CONCLUSIONS 

This  report  provides  a  numerical-analytic  solution  to  the  problem  of 
penetration  of  a  composite  skin  panel  with  a  slit  of  width  2b  by  a  quasi¬ 
static  magnetic  field  at  normal  incidence.  The  panel  is  modelled  by  an 
impedance  sheet  of  impedance  Zg.  The  basis  for  the  quasi-static  formulation 
used  is  provided  by  Latham  and  Lee  (1968).  The  solution  of  the  resulting 
dual  integral  equations  is  proqided  by  Sneddon  (1966).  For  this  particular 
problem  it  is  found  that  distribution  theory  and  the  generalized  Hankel  transform 
had  to  be  used  to  interpret  some  integrals  in  Sneddon’s  solution  procedure. 


A  convincing  argument  for  determining  the  values  of  two  constants,  A  and  C 
of  (8.7),  arising  in  the  formulation  needs  to  be  found.  The  values  for  these 
constants  have  been  chosen  so  as  to  allow  an  easy  solution  in  terms  of  Legendre 
polynomials.  C  has  been  chosen  to  be  one  of  two  different  values  depending 
on  whether  the  large  or  small  Zg  formulation  was  used.  For  the  special  case 
Zs=0,  an  analytic  solution  was  obtained  which  gives  the  same  boundary  values 
as  have  been  obtained  by  other  methods.  This  shows  that  the  basic  formulation 
is  correct. 

The  numerical  results  for  the  b  coefficients  of  (11.14)  and  (11.21) 

n 

were  disappointing.  It  had  been  hoped  that  b^  would  be  significantly  larger 

than  the  other  b  's  so  that  a  one  term  numerical  solution  could  be  found, 
n 

Unfortunately,  this  was  not  generally  the  case  and  problems  in  evaluating 
equation  (12.6)  for  n  >  1  did  not  allow  computation  of  the  potentials  in 
(10.9)  and  (10.10)  for  any  practical  cases. 

We  have  seen  that  the  Legendre  polynomial  solution  procedure  in  general 
did  not  adequately  model  the  physical  behavior  for  a  one  term  solution.  Perhaps 
another  solution  procedure  could  be  found  which  would  give  a  one  term  solution, 
but  It  is  not  clear  how  one  would  determine  such  a  solution  procedure. 
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